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The determination of the texture of the order parameter is important for understanding many 
experiments in superfiuid 3 He. In addition to reviewing the theory of textures in superfiuid 3 He-B 
we give several new results, in particular on the surface parameters in the Ginzburg-Landau region 
and bulk parameters at arbitrary temperature. Special attention is paid to separate the results that 
are valid at all temperatures from those which are limited to the Ginzburg-Landau region. We study 
the validity of a trivial strong-coupling model, where the energy gap of the weak-coupling theory is 
scaled by a temperature dependent factor. We compare the theory with several experiments. For 
some quantities the theory seems to work fine and we extract the dipole-dipole interaction parameter 
from the measurements. 

I. INTRODUCTION 

The superfiuid phases of liquid 3 He show complex behavior, which still can be understood theoretically. Many 
phenomena have been studied in a pure form in 3 He, and the knowledge can then be applied to other physical 
systems. For example, several structures of quantized vorticity have been seen in both A and B phases of 3 Hc [Q. 
The effect of impurities has recently been studied in many laboratories by aerogel immersed in liquid 3 He ||. Recent 
experiments on the Josephson effect show unexpected behavior || . Theoretical understanding of all these phenomena 
requires good quantitative understanding of the basic properties of superfiuid 3 He. This provides the motivation for 
the present paper. 

The purpose of hydrodynamics is to determine the behavior of a fluid on length and time scales that are long 
compared to some microscopic lengths and times Q. Hydrostatics is a subfield of hydrodynamics. It is restricted to 
study the equilibrium properties of the fluid. In simple fluids hydrostatics reduces to statements about the pressure 
variation in the fluid, which either rotates uniformly or is exposed to some external field. The problem becomes more 
difficult, if the fluid has some broken symmetry. Particular examples of these are liquid crystals and superfluids 4 He 
and 3 He. In both superfluids, the equilibrium mass current belongs to the scope of the hydrostatic theory. The order 
parameter of 3 He has also other degrees of freedom. The structure of those, which is often called texture, also has to be 
incorporated. In 3 He the hydrostatic theory is limited to length scales that are large in comparison to the superfiuid 
coherence length £o w 10 nm. The theory is valid at all temperatures. The hydrostatic theory can still be applied 
when the motion of the quasiparticles at low temperatures becomes ballistic rather than diffusive. The hydrostatic 
theory can be generalized to hydrodynamic theory by using conservation laws and adding the transport coefficients. 

Our purpose is to make a systematic presentation of the hydrostatic theory for the B phase of superfiuid 3 He. 
Part of the reason is that presently the various results are scattered over a large literature. Previous reviews treat 
hydrostatics only as a side topic and cover only a small part of the subject |5|-[ll). The existing papers are often 
unclear whether they treat general temperatures or are restricted to the neighborhood of the superfiuid transition 
temperature T c . We note that the hydrostatics of the A phase is better presented in the existing literature 11 



than the B phase considered here. In addition to reviewing, we present several results that have not been published 
before. 

We will begin with a general formulation of the hydrostatics of 3 He-B. Our approach is general enough to allow an 
external magnetic field and uniform rotation, both in the leading order. We write down an energy functional that 
consists of bulk terms and boundary conditions. All structures on the length scale of £o, such as surface layer or vortex 
lines, have to be treated as boundary conditions. The theory is found to split into two pieces: one for the superfiuid 
velocity and the other for the texture. The former is identical to the hydrostatics of superfiuid 4 He, whereas the latter 
can be solved only after the superfiuid velocity is determined. 

The coefficients of the hydrostatic energy can either be obtained experimentally or be calculated by some more 



microscopic theory. The calculation is discussed in sections [II and IV . The former considers the quasiclassical theory 
of 3 He [Q. The coefficients can be calculated using the weak-coupling quasiclassical theory at arbitrary temperature. 



* Associated internet page: http://boojum.liut.fi/research/theory/btex.html 
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We also discuss a "trivial strong-coupling" (TSC) model, where the weak- coupling coefficients are improved by scaling 
the energy gap. A different approach is studied in section |v|, where the hydrostatic coefficients are related to the 
parameters of the phenomenological Ginzburg-Landau (GL) theory at T w T c . 

Before any quantitative tests of the theory, we still have to determine the parameters that the quasiclassical theory 
needs as input. This is discussed in Section |y|, where we analyze several experiments. We find that certain quantities 
are reasonably well fitted using TSC model, but errors of 50% may occur for other quantities. 



II. HYDROSTATIC FREE ENERGY 



The superfluidity in a Fermi system arises from formation of Cooper pairs. A macroscopic number of pairs occupies 
the same pair state in the superfluid state [|l2| . The relative orbital wave function of a pair has p-wave symmetry in 
3 He and the spin state is a triplet f|[l3|]. The state of the pairs is thus described by an order parameter, which is 
a complex 3x3 matrix A a i. It gives the projections of the Cooper-pair wave function on the three p-wave orbitals 
(j>x, Pyi and p z , index i) and on the three spin triplet states (|— tt + II): i| IT + ll)i an( l I II + ll)i index a). In 
unperturbed B phase the order parameter has the form 

A aj = Ae^R aj (1) 

with real A, <f> and R a j. Here the amplitude A has a fixed (temperature and pressure dependent) value and R a i 
is constrained to be a rotation matrix, i.e. R a iR a j = 5ij. (Summation over repeated indices is assumed.) The 
phase <j) and the more detailed form of the spin-orbit rotation matrix R a i are not fixed on the scale of the superfluid 
condensation energy. These soft variables allow a dissipationless flow of both mass and spin. The order parameter 
can be interpreted as the wave-function of the center of mass of a pair. Using standard quantum mechanics, we can 
then define a mass-flow velocity 

v s = A V 0, (2) 
where m = 5.0097 ■ 10 -27 kg is the mass of a 3 He atom. In a similar way one can also define a spin-flow velocity 

Vs'a" — — ~a e a/37-R/3iVi? 7 i = Rai^ij /c-ffft/V Rf3k ; (3) 

4m 4m 

where Cijk is the maximally antisymmetric tensor. For example, if the axis of the spin-orbit rotation is constant and 
parallel to z, the pairs with spin states | ||), | ||), and | || + ||) flow with velocities v s + v^™, v s — VgP m , and v s , 
respectively. The three-dimensional rotation matrices are conveniently parametrized by an angle 9 and an axis n of 
rotation as 



Rij(n, 0) = cos OSij + (1 — cosd)nihj — s'm9eijkfik, (4) 

where n ■ n = 1. We note also the trivial identity a - R- a = cos 9 + (1 — cos#)(a • n) 2 , where a is an arbitrary unit 
vector. 

The soft variables <j) and R a i are determined by the interaction of the order parameter with various perturbations. 
The perturbations can be divided into external fields and boundary conditions. Experimentally, the most common 
field is the magnetic one H. It would also be straightforward to include the electric field, but we will neglect it here 
because its coupling is very weak jnj. The motion of the 3 He container can also be treated as an external field. In 
equilibrium the normal fluid component (velocity v n ) will follow the motion of the container, and the only allowed 
motions are uniform translation and rotation. The former is automatically taken into account because of Galilean 
invariance of the theory. The rotation will appear as a field V x v n , which equals twice the angular velocity. 

It is possible to construct a hydrostatic theory for any magnitude of the external fields. Here we assume the fields 
are small enough that the order parameter is not distorted strongly from the bulk form (Q). Strong perturbations, 
such as a surface or a vortex core, are treated as boundary conditions. These cause the order parameter to deviate 
substantially from the bulk form within a length scale of the coherence length £ ~ 10 nm, but on longer length scales 
also they act as weak perturbations. 

The degrees of freedom <j) and R a i differ crucially in the following respect. The mass flow can be written as 
j s = *p s ■ v s , where *p s is a phenomenological tensor. In the absence of external fields, *p s must be a scalar because of 
the isotropy of the unperturbed A ai (pi). The mass current has also to be conserved: V • j s = 0. Thus we arrive at 
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the Laplace equation V 2 </> = 0. Adding the boundary conditions, this completely determines 4>. The external fields 
cause only a small correction to this. Throughout the rest of this paper we neglect the small correction and start out 
from the assumption that the Laplace equation for <p is already solved, and thus v s is known. 

The problem that remains is to determine the rotation matrix R a i. What makes this different from <f> is that there 
exists interaction between the nuclear dipole moments of the 3 He atoms. It is of the form [jTJj 

F D . »,/^* + ««%) = 4A D y>rcos 9 (l + 2COS0). (5) 

Although this interaction is weak, it partly removes the degeneracy with respect to the rotation matrix. This means 
that the spin current is not conserved, but decays on a scale £d ~ 10 ^m. On the same scale, the rotation angle 9 
becomes fixed to arccos(— 1/4) sa 104°, which corresponds to the minimum of Fr> @, but the degeneracy with respect 
to the rotation axis n remains. Because no conservation laws exist, the rotation axis n is more susceptible to all kinds 
of perturbations than <fi. The subject of the rest of this paper is to study the texture, i.e. n(r) on a length scale ^> £d- 
We write down the free energy functional that governs the texture. The form of the energy terms is based on 
symmetry properties alone. The functional is valid in the limit of low fields and velocities, small gradients of the 
order parameter and weak coupling between the spin and orbital parts of the order parameter. The last condition is 
practically always satisfied because the coupling is due to the dipole-dipole interaction (^), which is small compared 
to the superfluid condensation energy by factor 10 _6 (1 — T/T c ) _1 . We neglect all constant terms, i.e. terms that do 
not depend on n. The leading terms in the expansion can be written as follows 
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F HV - -Ahv / d 3 r[H ■ R • (v s - v n )] 
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Fhvi = -Ahvi / d 3 rU ■ R • V x v n 
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The dipole-field term Fdh is discussed in Refs. the gradient term Fq in Refs. ||jl|[L7j , the dipole- velocity 

Fdv an d field- velocity Fhv in Ref- ||, and the first-order field- velocity term Fhvi in Ref. |lSf|~ Because of Galilean 
invariance only the combination v s — v n appears in Equations (0) and (^|). The superfluid velocity does not appear 
in the gyromagnetic term (0) because v s is curl free. Terms that are linear in n • H, n • (v s — v n ) or H ■ R ■ (v s — v n ) 
are prohibited by parity and time-reversal symmetry. Equations ((^)-(^l) serve as definitions of the parameters a, 
Adv, AhVj Ah~vi, Agi and Aq2, which depend on temperature and pressure. The calculation of these parameters are 
discussed in Sections HI and IV. For some of the parameters we use names given in Ref. fl7| instead of the more 
systematic names introduced here; for example Adh = a and Ash = d. 

The derivation of the gradient energy ( |To| ) deserves special consideration. Originally, one starts from a general 
expression that is quadratic in the spin velocity (|^). Making use of the properties of the rotation matrices, it is 
possible to simplify the energy to a form that is bilinear in the rotation matrices R a i- In addition to the two terms 
m ©, this for m contains a third term of the form diR a jdjR a i. The form ( p~C| ) can then be obtained by partial 
integration which converts the third term into the form diR a idjR a j . 

It should be noted that the partial integration of the gradient energy ( pL)| ) produces a surface term that is similar 
to Fqq below. Thus the value of the surface coefficient Age is unique only if the form of the bulk gradient energy is 
properly defined. Here the uniqueness of Age is guaranteed by restricting the bulk gradient energy to the form (|TJj|) . 

The gradient term can also be expressed explicitly as a function of n using the representation (0) |l7j . The needed 
identities are given in the Appendix. Our preference is to keep the shorter form (^) because a numerical algorithm 
can directly be based on it. 
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The dipole length £d is defined by £d = ^/Aq2/Ad- It is conventional to define dipole velocity vr> and dipole 
field -Hd by writing Ahv = 2a/(5up) and Adv = a ^D u rj 2 - We can also define a magnetic coherence length £h = 

65Ac2/(8ai? 2 ), which is inversely proportional to the field. The parameters defined here are temperature dependent. 
Near T c they reduce to constants that are commonly used. For example, £h — * R c Hb/H defined in Ref. fl7j| . 

In addition to the bulk terms d^)-jlO|), there are boundary terms. These energy terms originate from regions where 
the order parameter is strongly distorted from the form ([!]). We are here interested in two cases: surfaces and vortex 
cores. The boundary terms below are valid in the limit that the length scale of the distorted region (ss £o) is small 
compared to the dipole length £q. In reality this is well the case. It guarantees that the rotation angle 9 is not affected 
by the boundary. The form of the allowed boundary terms depends on the symmetry of the order parameter at the 
boundary. 

We assume that the surface structure has the maximal symmetry, i.e. time-inversion symmetry, rotation symmetry 
around the surface normal and reflection symmetry in planes perpendicular to the surface. (We note that also less 
symmetric states are possible pOfl.) We also assume that the curvature of the surface is small. Such a surface gives 
rise to the energy terms 



F S n = -d / d z r(H-R-sy (11) 
Js 

Fsuvi = -Ashvi / d 2 rll ■ R • s x (v. - v n ) (12) 
Js 

F SG = \ SG [ d 2 rs 3 R a ^ (13) 



s 



Or 



F SD = I d 2 r[6 4 (s-n) 4 -6 2 (s.n) 2 ]. (14) 

s 



Here s is a unit vector that is perpendicular to the surface and points towards the superfluid. The surface-field term 
Fsh is discussed in Ref. JlJ,[l7| , the surface-dipole term Fsb in Ref. Jl8|,|2l],|l7j , and the first-order surface-field- velocity 
term Fshvi in Ref. [fisf] . There are two contributions to the surface-gradient coefficient, Asg = ^sg + ^SG- The former 
comes from the equilibrium spin current that flows spontaneously along any surface p2f . In fact, the surface spin 
current = Xg G R a jeijkSk- (Note that J ss contains the factor h/2 for each fermion and thus has the unit J/m.) 
The other contribution Ag G comes from the partial integration that depends on the chosen form of Fq |i7f| . Note 



that there exists only one surface-gradient term (13) because R a iVR a j is antisymmetric in i and j. For the same 
reason the term ([l3]) does not depend on the normal derivative. We have constructed the definitions ©-(0) so that 
all surface (d, Ashvi, Asg, ^2, ^4) and bulk coefficients are non- negative, at least in the Ginzburg-Landau region. 

The order parameter is strongly distorted from the bulk form ([j]) in the cores of quantized vortex lines. Therefore 
the cores must be treated as boundary regions. We describe a vortex line by unit vector 1 that is parallel to the line 
and points in the direction of the circulation V x v s . The maximal point- symmetry operations of a vortex are a) 
reflection in plane perpendicular to 1, b) rotation around 1 (combined with a phase shift) and c) reflection in plane 
containing 1. The last one has to be combined with time inversion because otherwise the circulation would change 
direction. Assuming that the order parameter in the core has all these symmetries, we get the phenomenological 
terms 

Flu = Alh J d 3 r(H ■ i? • i) 2 (15) 
F LH i - A LH i / d'VH R \ (16) 



L 



Ald / d r [(I- n) + corrections]. (17) 



L 



The line-field term Flh is discussed in Ref. |23| ], the first-order line-field term Flhi in Ref- p4| , and line-dipole term 
-Fld in Ref. pa ]. Here L denotes the region where vortices are present. In _Fld only the dominant term is written 
explicitly. 
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It is well known that the vortex cores do not have the maximal symmetry E5|]. In the A-phase-core vortex the 
symmetry (a) is broken. Because this can take place in two different ways, we have to assign to each vortex line a 
new variable q that equals either +1 (left-handed vortex) or -1 (right handed vortex). This allows the line-gradient 
term M 

F LG = Alg J d 3 r ( q )i jRaj ^, (18) 

where (. . .) denotes the average because q may change from one vortex to another. Similar to the surface term Fqq, 
FjjQ arises from spontaneous spin currents. For an isolated vortex these currents form closed loops in the plane 
perpendicular to 1. All vortices in 3 He-B also have axial spin currents but they do not couple to external spin velocity 
in the lowest order because the net current vanishes. 

The double-core vortex also allows the term (|l8|). Additionally the circular symmetry (b) is broken leaving only 
discrete symmetry in rotations by tt. Thus an additional unit vector b perpendicular to 1 is needed to describe the 
vortex. This gives rise to line-anisotropy terms 

F LAH = A LAH ^rf 3 r((H-i?.b) 2 ) (19) 



■Flad = A LA d / d 3 r((h ■ n) 2 + corrections). (20) 



We note that there is flexibility in the definitions of the different terms. For example, the superflow around a vortex 
has to be counted into term Flh (^5|) in the region where the order parameter is strongly distorted but at larger 
distances it also can be included as the bulk term Fhv (||)- 



III. CONNECTION TO THE QUASICLASSICAL THEORY 

The energy terms (||)-(|2(]) contain a number of phenomenological coefficients. They should either be determined 
experimentally or calculated from a more microscopic theory than the hydrodynamic one. Pursuing the latter, there 
exists the quasiclassical theory Q. This theory bypasses the difficult many-body problem of strongly interacting 3 He 
atoms by concentrating in the low-energy range. It uses an expansion, where the relevant expansion parameter for 
the superfluid phases is the transition temperature divided by the Fermi temperature, T c /Tp ~ 0.001. The lowest 
nontrivial order in this expansion is known as the weak- coupling theory. It effectively contains the Bardeen-Cooper- 
Schrieffer theory as a special case, but it also reduces to the Landau Fermi-liquid theory in the normal state. This 
theory is adequate for some properties of superfluid 3 He, especially at low pressures, but it fails, for example, to 
stabilize the A phase. For many purposes it is important to continue the expansion to the next order in T c /Tp. We 
call this the strong- coupling theory. (Serene and Rainer use the name "weak-coupling plus", but we think this is too 
modest since there seems to be very little hope to calculate further orders in the expansion.) We will not go into the 
details of the quasiclassical theory, which is extensively discussed by Serene and Rainer jjj . 

It is important to realize that the quasiclassical theory is not microscopic in the sense that it would depend 
only on fundamental constants. Instead, it needs several parameters as input. This is especially a problem in the 
strong coupling case, which needs as input the scattering amplitude of quasiparticles (in the normal state) that is not 
accurately known. Additionally, the needed calculations are rather complicated at general temperature. There are two 
practical ways to proceed. The first is to restrict to the temperature region close to T c and use the Ginzburg-Landau 
theory. This approach will be described in Section |v| . The second way is to work at arbitrary temperature but to 
use the weak-coupling approximation in the quasiclassical theory. The latter approach is discussed in this section. 
At the end of this section we discuss how to improve the weak-coupling results by including a trivial strong-coupling 
correction. 

In the weak-coupling theory, the properties of the normal state are included via spin symmetric and antisymmetric 
Fermi-liquid parameters, F* and F z a (I = 0, 1, 2, .. .oo). We assume that the pairing interaction is effective in the 
p-wave channel only. The symmetry between particle and hole types of quasiparticles is consistently assumed in the 
quasiclassical theory. It turns out that all results presented below depend only on five Fermi-liquid parameters: F[ 
and F* with I = 0, 1, 2, and 3. (Infinite number of coefficients is needed in the hydrostatics of the A phase pq|.) In 
addition, the results depend on the mass density p of 3 He liquid, on the superfluid transition temperature T c and on 
the magnetic dipole-dipole interaction parameter go- 
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In the Bardeen-Cooper-Schrieffer model gu has the expression |HJ (in SI units) 



9 D = ^R 2 [h 1 N(0)nk B T 



(21) 



where R 2 is a renormalization constant and e c a high energy cut-off. (Note that our definition of g^> Q is different from 
that in Ref. [lq].) The Matsubara energies e n and the weak-coupling energy gap A(T) are defined in the Appendix, 
the gyromagnetic ratio of the 3 He nucleus 7 = —2.04 • 10 8 (T s)^ 1 , and the density of states at the Fermi energy 
2N(0) = (1 + |-Ff )(3m 2 p/7r 4 ?i 6 ) 1 / 3 . It is very convenient that the dependence of g^> on temperature is so weak that we 
can safely ignore it. In the weak-coupling approximation the constancy of gn{T) would be exact if the cut-off energy 
e c in ( pl|) were the same in the gap equation (|5l]) . (For the standard choice e c — ► 00 in the gap equation and R 2 = 1 , 
the relative variation of gn(T) is less than 10~ 5 .) In the trivial strong-coupling model (see below) Eq. ( fj"l|) gives the 
maximum variation at the melting pressure, where gn decreases monotonically by 1.3% when T decreases from T c to 
zero (assuming R 2 = 1). Because of uncertainties associated with e c and R 2 in Eq. (|2l|), we prefer to extract gr> from 
experiments, as will be discussed in section [v| 

For completeness, we give the results for nuclear magnetic susceptibility \ [PZI> superfluid density p s , and Ad (@) 



X = 2p N(0) 



+ a+lF$)Y 



1 + F a ( | + ly) + 1 F a ( 1 + (2 + F a )y) 



(22) 



1 - Y 



(23) 



Ar 



.9dA 2 



(24) 

1 - Z 3 (T) is the Yoshida 



All the following coefficients can be understood as corrections to these. Here Y(T) 
function, and the functions Zj(T) are defined in the Appendix. 

The basic principle for calculating the hydrostatic parameters is explained in Section VI of Ref. [Q. For the 
coefficient of the dipole-field energy -Fqh the main part of the work, the calculation of the gap distortion, is explained 
in detail in Ref. lEil. The result is 



5ffp 
2 



Foil 



I pa/1 
5^2 



F§)Y) 



5 - 



3Z B 



3F a Z 3 



^3 5(1 



(25) 



The coefficient of the dipole-velocity energy ([7]) can be calculated in a similar way and we obtain 



Adv = 5g D 



m vf 



F*Y 



1 - 



3Zs 
2Z, 



(26) 



Here m* is the effective mass given by m* jm = 1 + Ff/3. The Fermi velocity vf is related to basic parameters by 
vp = h(3Tr 2 p/m) 1 / 3 /m* . As far as we know, the expressions (25) and ( pq) have not been published before. Rather 
tedious calculation is needed for the coefficient in the field- velocity energy (||). This is done in Ref. |2j|, and we quote 
the result 



Ahv — 



m*/m 



A 2 (1 + ±FlYf 



±£7Mo(l + ^2 a ) 



1 + F$q + iy) + iF 2 a (i + (§ + f$)y) 



Zo 



9 r, 9 Z l 

— Z 5 H - 

10 10 z 3 



-z 7 



3F?Z, 



50(1 + \F-) 



(3Z 5 - 2Z 3 



(27) 



The gyromagnetic coefficient Ajjvi = because of particle-hole symmetry [g0[. The gradient energy coefficients are 
calculated in Ref. |3l]] , and can also be found in Appendix F of Ref. 0] . They are 



^G2 



h 2 p 



(l + ±Ff)(l + ±Ff)(l-n 



40mm* 1 + iFf(f + §Y) + iF 3 a (| + (f + \F?)Y) 



(28) 
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The structure of the order parameter near surfaces has been studied for a long time (for example in Refs. |32| ] 
and p2|), but the surface terms have been evaluated only quite recently. The gyromagnetic surface term ( p"2j ) 
vanishes identically because of particle-hole symmetry. As discussed above, the surface-gradient (|l3| ) term has two 
contributions: Asg = Ag G 4- Ag G . For the part that arises from the partial integration in the derivation of (|Io| ) 
we find Ag G = 2A G 2- The other part Ag G coming from spontaneous spin currents has recently been calculated in 
Ref. |33|. The same reference evaluates also the surface dipole coefficients in (|l4|). The field term ([ll]) has not been 
calculated. Until this is done, we can use an extrapolation of the Ginzburg-Landau result d = ^r(x n — x)£gl^0; where 
do = d/ 5hA 2 £gl is plotted by solid lines in Fig. |]. The Ginzburg-Landau coherence length £gl can be extrapolated 
to general temperature by £gl(T) = Hvp/"\/T0A(T). (Note that no strong-coupling correction to the weak-coupling 
A is allowed in this equation.) Xn is the susceptibility in the normal state [given by (|22|) with Y = 1]. 

For accurate calculation of the vortex terms, one needs a calculation of the order parameter in the vortex core. 
This has been done at general temperature by Fogelstrom and Kurkijarvi p4j| , but they do not give explicit values 
of Alh- However, at not too high rotation velocities, the most of the contribution to Alh comes from outside of the 
vortex core, and therefore a reasonable estimate at all temperatures is 

Alh » ^A HV (|v s - v n | 2 ) L « -^-fiA HV (ln — - 7), (30) 
2 2m tl 4 

where f2 is the angular velocity of rotation, R = yH/2mfi the unit cell radius of a vortex, and tl the radius of the 
vortex core. Because Alh is rather insensitive to tl, we may use r*L ~ £gl ||25[] , and use the same extrapolation of 
£gl as for the surface term. 

The weak coupling approximation used above is not expected to be accurate at high pressures, where strong 
coupling corrections are largest. As mentioned above, accurate strong-coupling calculations are very cumbersome, 
and introduce the scattering amplitude, which is poorly known. However, there is a simple procedure that is expected 
to take into account a major part of the strong coupling effects. This "trivial strong coupling correction" is to scale 
the energy gap A by a temperature and pressure dependent factor. This factor is tabulated by Serene and Rainer 
as a function of the temperature and the relative jump AC/C n of the specific heat at T c . The latter can be related 
to pressure according to the measurements by Greywall (B5[. Such scaling of A affects all hydrostatic coefficients 
©-© directly and/or via modification of the functions Zj and Y. 



IV. CONNECTION TO THE GINZBURG-LANDAU THEORY 



The hydrostatic theory was based on expansion of the free energy in small gradients and external fields. The 
Ginzburg-Landau (GL) theory is based on additional expansion in the amplitude of the order parameter [j36| . The 
expansion can be limited to a small number of terms near the superfluid transition temperature T c , where the order 
parameter is small. In most superconductors and in 3 He, the GL theory gives reliable results in the neighborhood of 
T c because the fluctuation range, where it becomes invalid, consists of a negligible temperature range just at T c . 

The order parameter in 3 He is 3 x 3 matrix A a j. The GL theory consists of writing down the terms in the free 
energy that are allowed by known symmetries. The superfluid condensation energy must be invariant in separate 
rotations in the spin and orbital spaces. This allows the leading terms [37[| 

-Fcond = J d 3 r[ — aA^A^i + PiA^A^A^jA^j + ^A^A^A^jAvj 

+fl3A* lli Al i A V jA tl j + P^A^AviA^jApj + p^A^A^A^jA^] . (31) 

The zero of the coefficient a = a'(l — T/T c ) defines the transition temperature T c . Other coefficients 0i (i = 1...5) as 
well as a' can be taken as constants in the expansion of -F C ond to order (1 — T/T c ) 2 . In the presence of nonuniform 
order parameter one needs the gradient energy f38fl 

2r j )(d i A fli )*d j A flj + [diA^fdiA^ + {2 n - 

l){d t A^)*d 3 A M + (cU w )*cU w - i(27? - l)e fcy -(V x v n ) fe ^A w -] (32) 



^G = K J d 3 r [(7 - 
= K J d 3 r [(7 
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with the Galilean-invariant derivative d = V + 2imv n /ft. In addition there are the energy caused by the magnetic 
field H |9|j|7| 



F H = J d 3 r (-ig m e Kl t V H K A*^A V i + g H H ^A^A^ + g' n H 2 A^A^) , (33) 

and the energy of the magnetic dipole-dipole interaction [fl5| 

F D =g D J cPriA*^ + A^A jt - fA*^). (34) 

We neglect all terms in the free energy that are independent of A a j . The gradient energy ( |32"| ) is parametrized using 
two dimensionless parameters 7 and 77, which are related to parameters introduced by Serene and Rainer |26j as 
7 = Kl/Kt and i] — Kq/Kt- The two different forms ( |32"| ) are equivalent, as can be verified by partial integration. 
In contrast to the hydrostatic case ([lO]), the surface term in the partial integration vanishes here because of the 
boundary condition SiA^i = 0. 

The parameters of the GL theory have been calculated using the weak-coupling quasiclassical theory, and the results 
are well known (see Refs. (i7j,|2j|, for example). There are two alternatives to incorporate the strong-coupling effects. 
One is to determine the coefficients purely experimentally. The /Vs, or at least some combinations of them, have been 
determined using experiments in the superfluid phases |4f|-^3]|. The other alternative is to consider the GL theory 
as a limiting case of the strong-coupling quasiclassical theory near T c |4^| . Here the problem of the poorly-known 
scattering amplitude is encountered again, but fortunately there exists model calculations for the most important 
coefficients. We give here a short summary of the results. 

There is a small correction to a' arising from finite lifetime of quasiparticles pq| . There are several suggestions 
for the Pi's that are based on different theoretical assumptions about the scattering amplitude and measurements 
in the normal state of 3 He [[f6|-[f8[. Although the strong-coupling corrections generally are small, they can be quite 
substantial in some combinations of /Vs. For example, /?345 = /?3 + /?4 + {3§ may differ 50% from its weak-coupling 
value. The corrections to K, 7, 77, gn, and g^ are calculated by Serene and Rainer in Ref. p6| . They find that 77 is 
unchanged from its weak-coupling value 1, but 7 increases from its weak coupling value 3 to w 3.1 at the melting 
pressure, g'^ is found to vanish even after strong-coupling corrections, and therefore it is dropped in the following. 
The parameter gm vanishes in the quasiclassical theory because of particle-hole symmetry, but this term is still kept 
because it is important in several situations. Its value is best extracted from measurements of the splitting of the A 



transition in magnetic field 49 50]. We have assumed that the nontrivial corrections to the dipole energy (|34|) are 



small, and therefore use the same coefficient gu as already discussed in Sect. III. 

The calculations in the GL theory are considerably simpler than in the general quasiclassical theory. Essentially all 
the hydrostatic parameters appearing in equations (^)-(|2^) have been calculated. We list below the bulk hydrostatic 
coefficients as functions of GL parameters. 



5ffDffH 

4/3; 



.545 



(35) 



5m 2 g D (7- 1)K 

A dv = ttt, I 36 ) 



•345 



Ahv = 2w2g £ 7 " 1)J " (37) 



h 2 f3; 



:i45 



, mgmh -4.r) + l)K , 

AHVl = fc , r, 77Z 7, — ^~ l J8 J 



h(p 4 ~ 3/3i - 35 ) 
aK 

2(3/3i2 + 0345} 



A G2 = KA 2 = ,, oa aK a , (39) 



T^=7-L (40) 

AG2 
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The first equality in (|39"| ) and ( (40| ) can be obtained trivially by substituting the B-phase order parameter dll) into 
the gradient energy The amplitude A of the order parameter is obtained by substitution into fcond ( ]3l| ) and 

minimization with respect to A. Equations (|3^)-(37) can be obtained by solving the GL equations in simple cases of 



axially distorted B phase. For example, the coefficient ( J35| ) can be obtained by first calculating the anisotropy of the 
gap due to a magnetic field and then evaluating the dipole energy for this gap. The gyromagnetic coefficient Ahvi 
( |38| ) has been calculated by Mineev (3(J. Because of deviation of 7 from 3, it is considerably larger than anticipated 
in Refs. 

Accurate determination of the surface terms requires a self-consistent solution of the order parameter near a wall. 
In the absence of fields, the order parameter A a i, which is normalized to unit matrix in the bulk, has real components 
A xx (x) and A yy (x) = A zz (x) near a surface located in the y — z plane. The surface coefficients are then obtained by 
integration: 

f°° dr ~ 

d = g H A^ GL (A 2 yy -A 2 xx ) (41) 

r°° flA 
Asg = KA 2 dx2( 1 -l)A vy ^ (42) 

f°° dx 5 

6 2 = . 9d A 2 Cgl/ — -(A 2 xx -6A xx A yy + 5A 2 yy ) (43) 
Jo Cgl 4 

f°° dx 25 

& 4 = 5dA 2 ^ gl / — — {A yy - A xx ) 2 , (44) 
Jo 4gl o 



where £gl = \J K/a. The surface term Fshvi can be found by calculating the order parameter in the presence of 
phase gradient: A a j(x,y) = Ai? Q ; exp(ifcy)Ay (2;). The coefficient is then given by 

Ashvi = ^shiA^l f %L \jmJ-[A xj A* yj - A yj A* xj \. (45) 

The surface coefficients d, 62, 64 [0, and Ashvi (ljj have been estimated before using simple models for the order 
parameter. We calculate them here by solving the order parameter numerically using the Sauls-Serene values f46|| for 
the coefficients /3j f20[| . For pressures below 1.2 MPa we smoothly interpolate the parameters to the weak-coupling 
values at zero pressure. We also assume the weak-coupling value 7 = 3. The calculations are done using boundary 
conditions appropriate for both specular and diffuse scattering of quasiparticles |j8) . (In the latter all the components 



of the order parameter have to vanish at the surface.) The integrals in Equations (41 )-(|§5j) (excluding the prcfactors 
(?h~A 2 £gl etc.) are dimensionless, and they are plotted in Figures [l] and ||. Because 62 > 264 > 0, Fqd ( ]l4| ) is 
minimized by n |j s [ pT| . 

All the vortex terms except Alg ( |l8| ) have been calculated in Ref. |25j . 

All the results of this and the previous section are, of course, identical in the limit where both theories are valid: 
weak coupling near T c . 

V. DETERMINATION OF PARAMETERS 

In his section we analyze a few experiments in order to deduce the values of parameters Fq, F|, and g D . We 



apply trivial strong-coupling corrections to the gap A, as explained at the end of Section III. For the molar volume 
v = tuNa/p as a function of pressure we assume the fit in Ref. p5| . 

The parameter Fq can be obtained from measurements in the normal state: the specific heat and the nuclear 
susceptibility. For the specific heat we use the measurements by Greywall |}5| . The susceptibility has been measured 
by Ramm et al J51J and by Hensley et al JE^] with essentially identical results. Unfortunately, it has been measured 
only below 2.9 MPa, and depending on the extrapolation alone, the relative error in 1 + Fg may be as large as 10% 
at the melting pressure. Examples are the simple fit Fq = —0.909 + 0.0055w cm~ 3 and the nonmonotonic Fq(j}) fit in 
Ref. [f33l. 
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12 3 

Pressure (MPa) 

FIG. 1. The integrals in Equations ©, (H), and (||). The solid lines give d/gHA 2 ^GL, dashed Xsg/KA 2 and dotted 
AsHvift/2m(?Hi A 2 £q L . Out of similar lines, the upper ones always correspond to specular surface scattering and the lower to 
diffuse scattering. 




2 3 

Pressure (MPa) 

FIG. 2. The integrals for surface-dipole terms in Equations ( [is) ) and (Q). The solid lines give &2/(?dA 2 £gl and the dashed 
&4/<?dA 2 £gl. For both quantities, the upper lines correspond to specular surface scattering and the lower to diffuse scattering. 
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The reduced nuclear susceptibility in the superfluid state, x(T)/Xn (|22|), has been measured by Corruccini and 
Osheroff , by Ahonen et al [[55| , and by Scholz et al [f56"| , [57| . (There has been a discrepancy between the susceptibility 
measured by NMR and by a SQUID magnetometer §13fl, but that is probably caused by difficulties in calibration p8|.) 
x(T)/xn depends only on two parameters, Fq and F-f. (We treat the trivial strong-coupling corrections as fixed, and 
use only the low-field limit of the Scholz data.) If we assume Fq given by Ramm and Hensley et al, only FJ 1 remains 
to be fitted. We find that a nonzero pressure-independent value of F% does not improve the fit essentially compared 
to Fg = 0. Since we believe that Ff cannot have strong pressure dependence, the simple choice F% = seems 
most attractive to us. Scholz finds FJ 1 ~ — 1 with a weak-coupling fit |57|], but this tendency is largely removed 
by the inclusion of trivial strong-coupling corrections. Note that the susceptibility data could be equally consistent 
with Ff = —0.7, say, but that would imply a systematic reduction of Fq by -0.025 from the results by Ramm and 
Hensley et al. Therefore we take Ff = in the following. For Fq we use the simple fit given above because it is in 
better agreement with the reduced susceptibility at the melting pressure Q than the fit by Halperin and Varoquaux. 
Note that at zero pressure our choice is not far off from the relation between _Fq _and F-f based on the " catastrophic 



relaxation" by Bunkov et al 59 1. There exists also other attempts to get F| 



The dipole constant gu has to be extracted from experiments because its value cannot be calculated accurately in 



the quasiclassical theory (Sect. III). The most straightforward way to get <?d is to measure the B phase longitudinal 
NMR frequency On . A direct measurement of 0|| has been made by Bloyet et al Q and Candela et al These 
experiments were done at low temperatures in the collisionlcss regime. According to the collisionless theory in a small 
magnetic field [Q 



h 2 N(0) \X 3 u 15 

where the function A(F) is defined in the Appendix. An alternative is to extract Qii from transverse NMR frequency 
in surface-oriented texture. These measurements have been done by Osheroff et al. |6q), by Ahonen, Krusius, and 



Paalanen |55[ (results tabulated in Ref. |66[), by Spencer, Alexander and Ihas 67 1, and by Hakonen et al. [p8|| . Because 
the external magnetic field in these experiments reduces the frequency difference between normal and superfluid 
precession, these experiments have to be analyzed using hydrodynamic theory. There On is related to g& via Jl5| 

ii x 



Both relations (46) and (|47j) imply that the temperature dependence of 0|| is fully determined by the energy gap A 
and the A function ([33|) or the susceptibility x (p2|). 

A third way to get <7d is so-called g shift of the transverse NMR frequency lo from the Larmor frequency uq. 
According to Ref. Q 

LU — LUq 4^qO 



ujo 5x 



(48) 



The g shift was measured by Osheroff |^] at the melting pressure and by Kycia et al |^] below 2.17 MPa. [We 
note that these measurements are done in such a high field that the expressions given in this paper are no more 
reliable near T c . However, Kycia et al measure the g-shift as a function of magnetization, and this plot is found field 
independent, both experimentally and theoretically |f42f . The only consequence is that the temperature T in Fig. |^ is 
the temperature according to the weak- field susceptibility (^2|), which differs from the true temperature near T c .] 

We plot <?d obtained by all three methods in Fig. ||. Let us first ignore the g-shift data (lines). It can be seen 
that the data for at each pressure is almost independent of temperature, as required by theory. We note that 
in order to reach this constancy it really is necessary that the longitudinal and transverse data of 0|| are analyzed 
with collisionless and hydrodynamic theories, respectively. Equally important is that we use trivial strong-coupling 
corrections. The value of grj (but not the temperature dependence) also depends on Ff and T c , for which we use the 
measurements by Greywall pq] . 

The go data as a function of pressure is plotted in Fig. || It contains all the same data as Fig. || and some additional 
data. It can be seen that the longitudinal and transverse measurements of On agree very well at intermediate pressures, 
but there is a difference at both high and low pressures. Theoretically the ratio of collisionless and hydrodynamic On 
(at the same pressure) depends only on the reduced temperatures T\/T c and TijT c of the two measurements and on 
Fq, F£, and AC/C n . The difference in the effective <?d obtained by the two methods seems at low pressures much 
larger than the expected uncertainties of the parameters, and remains unexplained. Both the positive slope of <?d(F) 
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0.0 



_1_ 



0.0 0.2 0.4 0.6 0.8 1.0 T/T Q 

FIG. 3. Plotted as a function of temperature, go interpreted from different experiments according to trivial strong-coupling 
model. In ideal case go should be independent of temperature. We have used the transverse NMR frequency measured by 
Osheroff et al. Q (A 3.44 MPa), by Ahonen, Krusius, and Paalanen || (0 3.2 MPa, T 2.9 MPa, V 2.54 MPa, ■ 2.11 MPa 
+ 1.87 MPa,), and by Hakonen et al. || (• 2.5 MPa, o 1.55 MPa, ♦ 1.02 MPa, x 0.5 MPa, A 0.05 MPa). We also used the 
longitudinal NMR frequency measured by Candela et al. |^l[ (□, from top to bottom 3.3, 2.1, 1.2, 0.61, and 0.03 MPa). For 
comparison, there are curves that are determined from the g-shift measurements by Osheroff [^s| (curve with label x5, reduced 
by factor 1/5) and by Kycia et al Bl (curves, from top to bottom 2.17, 1.3, 0.7, 0.3, and 0.11 MPa) 




12 3 

Pressure (MPa) 

FIG. 4. Plotted as a function of pressure, go interpreted from different experiments according to the trivial strong-coupling 
model. The data from f2|| in transverse NMR (+), from longitudinal NMR by Candela et al (□), and g shift (bars) are the same 
as in Fig. ^[ In addition there is data from zero-field longitudinal NMR measurements by Bloyet et al (0) and transverse 
surface-oriented NMR measurements at saturated vapor pressure in the temperature range T/T c = 0.83 . . . 0.86 by Spencer et 
al |j7| (x). The solid and dashed lines approximate the data from longitudinal NMR and g shift, receptively. The dotted line 
is the model of Eq. (^) with e c = 0.45/cbTf. It is argued in the text that go given by Qy is the correct one, and the difference 
between <;d's deduced from the g shift and flu is a measure of nontrivial strong-coupling corrections. 
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(Fig. ||) at high pressures and the difference between transverse and longitudinal data could be reduced by giving Ff 
a negative value (« —0.7 at high pressures), but only at the expense of impaired fit of the x(T)/Xn- 

We believe that the nontrivial strong coupling corrections are small in the expressions for f2y (^6|)-(p7|), which are 
based on expectation values in an unperturbed order parameter. [There are nontrivial corrections to x P^l i but as 
long as (EEj) can reproduce (possibly with incorrect FJ 1 ) the measured x, the value obtained for #d is unaffected.] 
The collisionless f2|| data is fitted by solid line in Fig. |]. The values obtained for gn depend on Fq , T c , iV(0), 
F%, and AC/C n , and must be revised if more accurate values of these become available, for example, via improved 
measurement of the temperature |70| . 

We can also estimate <?d based on the simple model of Eq. (|2l]) . Assuming R 2 — 1 and making the sum at T = T c 
gives 

2D = Sf^V(0)lnH^y. (49) 



40 V k B T c 

We take the cut-off energy e c proportional to the Fermi temperature defined by Tp = 3p / '4N (0)&bTO. As shown by 
dotted line in Fig. ^, the resulting expression fits nicely the experimental data for the constant of proportionality 
Ec/^bTf = 0.45. The agreement may be accidental, however, because there is no fundamental justification for the 
approximations made. 

Let us next study the g^ data based on the g-shift. It is also rather independent of the temperature (lines in Fig. |^). 
Fig. |^ shows that the gjy data (bars) at different pressures are well consistent with each other: the low pressure data 
extrapolates well to the melting pressure data by Osheroff. However, g-Q deduced from the g shift differs essentially 
from the determinations based on fly, especially at high pressures. The reason for this is that the expression for a 
( p5j ) has substantial strong-coupling corrections that are not included in the scaling of the energy gap A. This can be 
seen by comparing the T — > T c limit of trivial strong-coupling a (|25|), denoted by axsc, with the Ginzburg-Landau 
limit a GL (Ph. We find 



owe 

QGL AGL ffH P345 , Kf y> 

~ 3vc « — ■ ( 50 ) 



otsc awe 5h # 



•:t 15 



where WC denotes weak-coupling. It is well known that ^345 differs substantially from its weak coupling value 
|j46|,ffl 43 1 . This explains the difference in the apparent gu deduced from g shift and fiy. The new thing in the present 



analysis compared to Ref. J42J is that the difference is not limited to the Ginzburg-Landau region, but because of the 
weak dependency of the apparent gr> on temperature (Fig. ||), it persists almost unchanged at all temperatures. 

We conclude this section with a comparison of theoretical and experimental dipole velocity ud- Theoretically this 
quantity is related to coefficients a ( p5| ) and Ahv ( p7| ) by the relation v 2 ^ = 2a/ (5 Ahv)- It has been measured by 
Nummila et al J7l| . Originally they compared their result to a theory that turned out to be in error, see discussion in 
Refs. |2^,|7|]. The comparison with the present theory is given in Fig. We have used trivial strong-coupling theory, 
parameters as described above and <7d from solid line in Fig. m 



Because both a (35) and Ahv ( p7| ) are proportional to P^ i5 in the Ginzburg-Landau region, the uncertainty discussed 
in connection with a is expected to cancel out in vu- There also exists a direct measurement of Ahv p9| - It also 
shows deviation from the trivial strong-coupling model, but the differences are not of similar type as for a, and are 
presently not understood. 

Above we have discussed all input parameters of the trivial-strong-coupling hydrostatic theory except Ff and F£ . 
Out of the bulk terms only the gradient coefficients ( p8| ) and ( p9j ) depend on these. There is several independent 
evidence that F* w — 1 at high pressures ]73|-|75|] but F£ is not known. 



VI. CONCLUSIONS 



We have presented a summary of the hydrostatic theory in superfluid 3 He-B. Several new analytic and numerical 
results were included. Some experimental data was analyzed in order to extract the parameters of the theory. A 
particular goal was to understand how well the B phase is described by the trivial strong-coupling model. We found 
that some quantities (x, <?d, ^d) can successfully be calculated, but there are other quantities (a, Ahv) that may be 
wrong by 50% in this model. The parameter g^> has direct relevance also for the A phase, where it has been used in 
comparison between theory and experiment (see Ref. [79] , for example) . 

The first application of the results calculated here has been the comparison of the ratio d/a [from equations ( |4l| ) 
and d35|)] to experiments in Ref. f?% . More recently, Kopu et al [M applied the hydrostatic theory to a rotating 
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cylindrical container. They calculated the NMR line shape and studied the optimal conditions for observing single 
vortex lines. Another application is the Josephson it state observed recently ||. This effect was found to depend 
crucially on the texture at the Josephson junction j^,^. The texture is also essential in several experiments of 
superfluid 3 He in aerogel. For example, the identification of the B phase was based on its texture-dependent NMR 
spectrum p3]. For the present, the textural parameters in aerogel have been evaluated only in the Ginzburg-Landau 
region in the homogeneous scattering model Q] . These developments demonstrate that there still are open problems 
in superfluid 3 He and in many cases a proper understanding of the hydrostatic theory is a prerequisite for solving 
them. 
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APPENDIX 



We give here some equations that complete the theory presented above. The weak-coupling energy gap A(T) is 
determined by the equation 



In 



T 

Tc. 



71— — OC 



1 



1 



= 0, 



where the Matsubara energies e n = 7rT(2n — 1) with n = 0, ±1, ... ± oo. The Zj(T) functions are defined by 
Y(T) = 1 - Z 3 (T). The A(T) function J 

oo 



which equals to 1 — f(T) defined in Ref. j82|], can be written as 
A 



74 + tf(^+tf + A) 

The numerical calculation of the functions is discussed in Ref. [p3| . 

The gradient energies @ and @ can be written in different forms using the identities jjj 



(51) 



(52) 



(53) 



(54) 



diR aj diR aj = 4(1 - cos9)(d l h j ) 2 = 4(1 - cos(9){(V x n) 2 + (V • n) 2 + V • [(n • V)n - n(V • n)]} 
d,R ai djR aj = (1- cos6»)[2(V x n) 2 + (1 -cos0)(V • n) 2 - (1 -cos0)(n- V x n) 2 - 2sin6>(V • n)(n- V x n)] (55) 
d t R a:j d 3 R ai = (1 - cos6»){2(V x n) 2 + (1 - cos6»)(V • n) 2 - (1 - cos6»)(n • V x n) 2 

-2sin0(V • n)(n • V x n) + 2V • [(n ■ V)n - n(V ■ n)]} (56) 
SiRaidjRaj = -(1- cos0)s- [(n • V)n - n(V • n)]}. (57) 

Using these it can be seen that Fq ( |l0| ) has a pure divergence term cx V • [(n • V)n — n(V • n)]. The prefactor of this 
term is half of the value of that by Smith, Brinkman and Engelsberg [0. With present definitions the other half is 
transferred to the surface term 
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FIG. 5. The dipolar velocity «d- The circles are experimental data from Ref. |7l{ ]. The solid lines are theory based on 
Equations and (VA), where all parameter values are fixed by other measurements as explained in the text. 
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